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1 Introduction 

Theoretical study of the reaction p+p — > j + c+X 
is interesting for several reasons. The two particle 
inclusive cross section probes QCD dynamics in 
more detail than a single particle inclusive cross 
section. Specification of the cross section in next- 
to-leading order in QCD is correspondingly more 
challenging, and data will permit a study of var- 
ious correlations including those in the rapidities 
and transverse momenta of the final photon and 
charm quark. Second, the data may permit a di- 
rect measurement of the charm quark density of 
the proton if it can be shown, as we do Em, that 
the hard scattering cross section is dominated by 
the contribution from the eg scattering subpro- 
cess, where the incident c quark is a constituent 
of one of the initial hadrons. Third, it is valuable 
to be able to compare data obtained by the CDF 
Collaboration U on prompt photon production in 
association with charm decay products, such as 
e^- , fi^ and D^, with next-to- leading order pre- 
dictions. | 

We have published two next-to-leading order 
calculations in perturbative QCD ofr-the reaction 
p+p — > 'y + c+X at collider energiesSla. In most of 
the data of interest, the charm quark carries a rel- 
atively large value of transverse momentum, p^, > 
10 GeV. Since this value is large when compared 
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Figure 1: Cross section as a function of z at %p = in 
leading order. 



to the mass of the charm quark, we neglect the 
mass of the quark, and we include hard-scattering 
subprocesses in our calculation in which the charm 
quark is a constituent of the initial hadrons. For 
values of p^, less than 10 GeV or so, it would be 
necessary to retain a finite mass for the quark and 
to examine a different set of subprocesses in which 
the massive quark is not a constituent in the ini- 
tial state but is, instead, produced in the hard 
scattering □. A calculation with m c ^ also re- 
quires resummation of logarithmic terms of the 
form a™ln"(py/m c ) for large values of px- 

In leading-order, the contributions to associ- 



ated production include the direct term c + g — > 
7 + c along with several photon fragmentation pro- 
cesses typified by the hard scattering subprocess 
c + q — > c + g followed by long-distance fragmen- 
tation of the final q into a photon: q — > 7 + X. 
In next-to-leading order, there are gluon-exchange 
loop corrections to the lowest-order direct process, 
real gluon emission radiative corrections to the 
lowest-order direct process, and several new three- 
body final-state contributions: g + g^c + c + j, 
q+q — ► c+c+7, c+q — > c+g+7, c+q — ► c+g+7, 
c + c — > c + c + 7, an d c + c — > c + c + 7. We include 
all of these. In principle, we should also include 
fragmentation to next-to-leading order. However, 
in view of the small size of the fragmentation coru 
tribution and the theoretical problems identified cl 
with isolated prompt photon production in ncxt- 
to-leading order, we treat the fragmentation terms 
at lowest order. 

Xhe distinction between our two calcula- 
tions M resides in the different techniques used 
in performing the phase-space integrals. In the 
first paper, we insist on going as far as we can us- 
ing purely analytical techniques. In the second, 
we use a combination of analytical and Monte 
Carlo techniques 0. The second approach allows 
implementation of isolation cuts and other exper- 
imentally relevant selections. In both cases, we 
evaluate the gluon exchange and the real gluon 
emission diagrams analytically, expose the infra- 
red and collinear singularities, cancel the infra-red 
singularities, and factor the collinear singularities 
into incident parton densities or final-state frag- 
mentation functions. We use dimensional regular- 
ization to deal with the singularities. Although 
somewhat restrictive, the purely analytic calcula- 
tion EJ provides the size of the two-particle inclu- 
sive cross section, the relative importance of var- 
ious contributing subprocesses, and an important 
check on the more versatile combination □ of ana- 
lytic and Monte Carlo methods. 

Some of our analysis is done in terms of the 
variable z — -p^ ■ Pt/(Pt) 2 ' ano - some of the 
results are presented as histograms in z for fi- 
nite intervals in z. If z is positive, the photon 
and charm quark are in opposite hemispheres. If 
z = 1, the photon and charm quark have equal and 
opposite values of transverse momentum. This 
situation occurs for the leading-order direct pro- 
cess, for gluon exchange loop corrections to the 
leading-order direct process, whenever an emitted 
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Figure 2: Cross section as a function of z in next-to-leading 
order. The dominant contributions are shown. 



gluon becomes soft, and when a third parton in a 
threc-parton final-state process becomes collinear 
to an incident parton. The point z = 1 is associ- 
ated with various soft and collinear poles, exposed 
when the cross section is expanded as a function 
of z. The singularities are expressed in terms of 
delta functions, 5(1 — z), and "plus" -distributions, 
1/(1 — z) + . The singularities are integrable, and 
results are expressed as integrals over finite inter- 
vals in z. 

2 Numerical Results 

Our results are presented at a center-of-mass en- 
ergy y/s = 1.8 TeV. The renormalization and fac- 
torization scales are taken as /1 = pip. We sum over 
charm and anticharm production throughout. In 
Fig. 1 the net lowest-order contribution is shown 
as a function of z at p T — 15 GeV. The lowest 
order cross section is made up of the lowest order 
direct term eg — > 7c, proportional to 5(1 — z), and 
fragmentation contributions. In Fig. 1, the direct 
term provides the peak at z = 1, and the various 
photon fragmentation terms contribute in the re- 
gion z > 1. Note the absence of any contribution 
in the region z < 1, an unrealistic prediction that 
illustrates the inadequacy of the lowest order ap- 
proach. The cross section in z resulting from our 
the next-to-leading order massless-charm calcula- 
tion is displayed in Fig. 2. The next-to-leading or- 
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Figure 3: The transverse momentum dependence of 
dcr I 'dpTpdy* 1 dz , for z integrated over the interval 0.2 < z < 
2.0. The upper solid line shows the sum of all subprocesses 
through next-to-leading order. The dashed line shows the 
sum of the 0(a a ) and O(a^) contributions from the eg 
initial state. The O(a^) contributions from the gg and cq 
initial states are shown as dash-dot and dotted curves. The 
lower solid line shows the O(a^) contribution from the qq, 
cc, and cc initial states. 



Figure 4: Cross section dcr/dp^j, as a function of p^, with 
15 < p7p < 45 GeV, and the rapidities of the photon and 
charm quark restricted to —1.0 < y < 1.0. The solid curve 
shows the cross section with no further restrictions and 
the dashed curve the result after the additional selection 
z > 0.1. The dotted and dash-dot curves display the results 
after photon isolation restrictions are applied, in addition 
to the cut on z. 



der collinear initial-state contributions reduce the 
peak at z = 1, and the next-to- leading order three- 
parton final-state contributions broaden the distri- 
bution. The cross section is finite at all values of z, 
similar to the situation observed in experiments. 
In Fig. 2 we also present the individual contribu- 
tions from the most important subprocesses. The 
eg initiated process dominates the cross section, 
but there are important contributions from the gg 
and cq initiated subprocesses in the region of small 
Pj,. The predicted distribution in pip is shown in 
Fig. 3. As might be expected, the contribution 
from the qq initial state begins to become impor- 
tant at large pip. 

Distributions in p^ and in the rapidities of 
the 7 and c quark may be found in our publi- 
cations Era. In the second paper, we also discuss 
the influence of photon isolation, and we compare 
the magnitude and shape of the next-to-leading 
order and leading order cross sections for various 
observables. As an example, we present Fig. 4, a 
distribution in p^ for different selections on z and 
the isolation cone variable R. 

Positive correlations in rapidity at collider 



energies are a dynamical property of the hard- 
scattering matrix element Q. Indeed, it may be 
shown Q that the parton-level matrix element is 
proportional to cosh -1 (Ay), where Ay is the dif- 
ference in the rapidities of two final-state particles 
in a two-particle inclusive reaction. In Fig. 5, we 
present the distribution in the difference of the ra- 
pidities of the photon and charm quark. In this 
example, the 7 is selected to have rapidity in the 
forward region 1 < y 7 < 2. Under these circum- 
stances, as may be observed in Fig. 5, the charm 
quark is predicted to have positive rapidity also, 
with Ay = y 7 — y c peaking a value only slightly 
larger than 0. 

In Fig. 6 we compare our results to CDF datai 
for photon plus ^ production, where the (ji^ is 
identified as the decay product of a charm quark. 
The three upper points are obtained 19 from the 
Monte Carlo event-generator Pythia whereas the 
lower point is our theoretical expectation. The 
Pythia cross section lies substantially below the 
data while our cross section exceeds the data some- 
what. 
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Figure 5: Cross section dcr/dAy as a function of the dif- 
ference Ay = y 1 — y c of the rapidities of the photon and 
charm quark, for p + j)-*7 + c + Jfat ^fs = 1.8 TeV. 
The ratio z is restricted to z > 0.1, and photon isolation 
is imposed, with R = 0.7. The transverse momentum and 
rapidity of the photon are selected to be in the intervals 
15 < < 45 GeV and 1.0 < y 1 < 2.0. The dashed curve 
shows the behavior at leading order. 



3 Conclusions 

We presented the results of next-to-leading or- 
der QCD calculations of the inclusive production 
of a prompt photon in association with a charm 
quark. The two-particle inclusive cross section al- 
lows more refined tests of the underlying QCD 
dynamics, revealed, e.g., through rapidity corre- 
lations. Our cross section shows the expected 
reduced dependence on the choice of the fac- 
torization/renormalization/fragmentation scale in 
QCD. The sensitivity of the cross section to the 
charm quark content of the proton is preserved 
in next-to-leading order. A comparison with the 
CDF data shows reasonable agreement. Further 
work is underway to include effects associated with 
the finite mass of the charm quark and with charm 
fragmentation: c — > DX or c ^ /iA. In other 
work, we are investigating associated production 
of photons and charm quarks at HERA energies 
and the associated production of photons and bot- 
tom quarks at collider energies. 



A Pylhia (LO+Paiton Shower) 
• NLO QCD 



T 

9- 15 

T 



ye 



ye YH YD 

yc Cross Section Measurement Type 



Figure 6: Ratio of the measured ■ycX cross section to that 
predicted, for various charm decay channels. 
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